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Abstract 

The effect of resonant tunneling on magnetoresistance (MR) is studied theo- 
retically in a double junction system. We have found that the ratio of the MR 
of the resonant peak current is reduced more than that of the single junction, 
whereas that of the valley current is enhanced depending on the change of 
the discrete energy-level under the change of magnetic field. We also found 
that the peak current-valley current (PV) ratio decreases when the junction 
conductance increases. 

I. INTRODUCTION 

As nano-fabrication technology using magnetic materials advances, the magne- 
toresistance (MR) in mesoscopic systems has attracted growing interest, mainly be- 
cause of its possible applications in storage devices The resonant tunneling 
phenomenon is one of the events expected to occur in magnetic nano-particles or 
thin film systems where discrete energy-level structures are prominent. Ferromag- 
net /insulator /ferromagnet /insulator /ferromagnet (FM/I/FM/I/FM) is a basic structure of 
this system. The electron tunneling through the insulator is considered to be based on the 
independent polarized electrons and the difference in MR is considered to be derived from 
the difference of the spin polarized free electron tunneling Recently, Zhang et al. |^ 
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calculated the I-V characteristics of the FM/I/FM/I/FM structure based on the Tsu-Esaki 
formulation (S-matrix theory) and showed that the MR is more than 90%, which is a great 
enhancement compared with the case of a single junction . However, the origin of the en- 
hancement is not clear. In addition, although the peak current-valley current (PV) ratio is 
one of the key factors of resonant tunneling, the effect of the MR on this ratio has not been 
clarified. This is because S-matrix theory is not easily applied to analyze the underlying 
physics, in spite of its usefulness as a more realistic method for self-consistent calculations. 
In this paper we discuss the effects of the resonant tunneling on the MR analytically, based 
on the two band spin polarized free electron model. We assume that the capacitances of 
two junctions are large and neglect charging effects. The general formula of current-voltage 
characteristics is obtained by elaborating that derived by Jauho et al. and is shown to be 
a useful formula for investigating the detailed physics of a double barrier structure. As we 
are using Green's function method, our model can be extended, by refining the self-energy 
part of Green functions, to be applied to the case in which free electron approximation can- 
not be used because of the existence of some scatterings, although it is not the subject of 
this paper. 

II. FORMULATION OF THE DESCRIPTION OF THE SYSTEM 

The Hamiltonian consists of the electronic part H-^, and the transfer part H^^. The 
electronic part consists of electrode and island parts, 

ao^~^ ^~^ EmgdlnfT^ma 1 (l) 
k,aGL,-R,(T ma 

where a represents a set of parameters that, together with wave vectors k, completely specify 
the electronic state of the left (L) or right (R) electrode, and m specifies the energy levels of 
the central island. By an internal magnetic field and with the Pauli spin matrix a, the 
energy dispersion relation is expressed as 

E^aa = h^K/{2m) -K-a. (2) 
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Hereafter we write Ei^^ai = /?-^k|/(2m) — and -Ekaj, = ^^k^/(2m) + ha- The transfer part 
is described by 

-^T = [Vnkca{t)cl,aadna + h..C.]. (3) 

The current at junction a (=L or R) is given by 

— ^~^^Ji IVrikaa (Ckn-rr'^rag-) ~ Kikao-('^ri(TCkao-)] 
nkcr 

= (-f|Re|5:Kk..(t)G4„,(t,t)|, (4) 

where /3=0 for the left junction and P=l for the right junction and G',^,^^^(t, t') = 
'^{^i.aai'^)'^na(t')) is an analytic continuation of the contour-ordered Green's function 
G^j^Q,^(r, r') which is defined in the interaction representation by 



xexpi^--J^HT{ri)dnj}), (5) 

where Tc is the contour-ordering operator. We assume that electrons in the left and right 
electrodes are noninteracting. Then the only nonvanishing terms in Eq. are those in 
which cj^Q(T') is contracted with Cka(Ti) in the exponential term. We then obtain, after 
analytic continuation to real time, 

+G<_(t,ti)^7L.(ii,0], (6) 
where G^^^{t,t') = 'i{dlnait')dna{t)) is the Green's function of the central island, and 

9LAti,t2) =i{claAt2)cka{ti)), (7) 
9LMut2) = -^{c^aAtl)cLAh))■ (8) 

Substituting Eq. (||) into Eq. @), we obtain 



Ut)= (-1)^-2 E / dt'V:,^^„{t'){Gl^^M 



where G^,^^{t,t') = ~i{dnait)dl^^(t')) . By introducing a noninteracting self-energy, Sg<, 
defined by Gq^ = GqSq ^Gg, and using the Dyson equations, (1 + G^T,'^)Gq = G"" and 
^^(1 + S^G") = G'^ we may cast the Keldysh equation G>< = (1 + G'1:')G><{1 + S^G") + 
Qrj^Qa jj^^Q following fomi [^]: 



xG::^„.(t2,o, (10) 



where 



^to^nin.2a('^l)^2) — SQ^^„2a('^l)^2) + Sx'iinj^ (^1 , ^2) (H) 

The first term on the right-hand side describes the self-energy of the island in the absence of 
disorder, interaction, and tunneling. Since it corresponds to the free part, it is infinitesimal: 
^?[/o(e) - 1/2 =F 1/2], where r] is a positive infinitesimal, > (<) refers 
to the minus (plus) sign, /o(e) = (exp (3{e — Vd) + 1)~\ and Vd is the bottom of the island 
energy. The second term describes a self-energy due to tunneling: 

^Tnin2(Tln) '^2j - 2^ p S'koa l^^l ) ^2 J • 

ko 

Because of the infinitesimal factor r], the free part is important only when the remaining 
parts are absent (it is not the present case). 

In the following we focus on the case in which r„„/kao-(t, t') = 2ti V*]^^^{t) Vn'\^aa{i') 
is a real function of t - t' . From Eq. (D, J^-Jr =eE.IZodt' {G>^^{t,t')E<^^^{t' ,t) - 
GnmAt,t')^Tmnait',t))- Bccausc G><^{t,t') coutalus S^<„^(t,t') (Eq. (|T2|)), the conserva- 
tion of current through the two junctions is automatically satisfied, that is, Jl = Jr- Thus 
we can discuss the current through the junctions by Jl- 

The retarded and advanced Fourier-transformed Green's functions at the central island, 
Gnn'ai^) G'^^,^{e), are derived from the Dyson equation: 



= e-{Ena + Vda)-Kn'a{e)±^ (2r/+r„„v(e)) , (13) 

where r„„/<^(e) = 2ImE:^.„„,^(e), and 27r/iST^„,^(e)=Eka r„„'kQ<7(e)5'ka<T(e)- The real part of 
the self-energy A„„/o.(e) shifts energy-levels in the central island and we will regard this effect 
as included in our assumed one-body energy levels of the quantum dot. Here r„„/o.(e) shows 
the half-width of resonant peaks : 

rnn'a(e) =Jlj d(^l'^nn'ka{(^l)S{e - €1 - Ekaa) (14) 
ka 

With A^^n2aie) = ^(G;^n2^(e) -^^^^.^(e)), and henmA^) = G;„^^(e)G'^^^^(e)/A„,„2<7(e), Jl 
can be cast into the following form, 



kk'nm 

nin2<T 

X Qnma (c) {/l (^kLa) " /r (^k'/Ja) }^ (c " E^La) ^ (e - ^k'i?^) • 



(15) 



This is a general expression of current where /L(e)=l/(e'^*^'^~^^'^~^^^l) and /i?(e)=l/(e'^(^~^^''^ 
1) with Epa being the Fermi energy of the electrodes. Hereafter, we apply this expression 
to simple cases. First we assume that S^'^^(e) = dmn^n"""- This corresponds to a situation 
in which energy levels in the island are mutually uncorrelated during the tunneling process. 
Then all Green's functions are diagonalized and ^nin2o-(e) reduces to 

^"'^'^ " [e - {E^, + V^) - K^e)]' + [T^Ae)f /4" ^^^^ 
and ©„o-(c) reduces to r„o.(e)~^. 

rnn'a-(e) IS related with the density of state (DOS) of the electrode such that 

nn'aa (e) - ^nkaa^n'kaaS{(^ " -^kaa) 
k 

-^^^K:k^aaK'k^aa5(e - ^kaa) 
= KkFaa^WkFaaDaai^), (17) 
5 



where Daa{E) is a density of state expressed by Daa{E) = ^ (f?) ^ VE + aha in a three- 
dimensional space. This approximation will be most suitable when the interfacial material 
is a quantum box or particles. From the relation Taa{^) = 27rDQ,o-(e)|VkcKTp, resistances 
i?ao-(tt = L,R) can be evaluated to be Raa/Rx = (Xaa{Epa-)Dd{Ep^))~^, where Rk is the 
resistance quantum h/e"^ =25.8kr2 and Dd{Epcr) is the DOS of the island at the Fermi energy. 
In this paper MR is discussed in terms of the change of Taa instead of Raa 

We see r„„o-(-EkQcr) as a function of the internal magnetic field : r„Q,o.(-EkcKT) = ^naaiha). 
Thus we have the resonant tunneling formula under an external magnetic field, H, as, 

Jl{H) = fdE]^La "^'"^ "'^p^ .^j. 7— Tal/il^kLa) - /ij(^kLa)}, 

(18) 

where Sq = ma.x{eV — ah^, -fr/i/j). 

Here we compare the MR of the resonant tunneling current with that of the single 
junction by their conductances. The conductance of the system is given by = dJ^/dV, 
and the magnetoresistance is defined as 

Gil - G^i 



MR^ 



We have a picture that the change of magnetic field makes the distribution of the DOS 
(Fig.0) and the Taa{ha) (Eq.(p^). We consider the derivation of Taa{ha) from the case 
where the internal magnetic field is zero, = 0: 

TnLaihL) = r„L.(0)(l + Ai,(/lL)), 
^nRa (M =rn/?a(0)(l + A^,,(/lfi)), 

= K.(0)(1 + 7.(/in)), (20) 

where \ Aaaiha)\ ^ 1 {a=L, R) and |7cr(/in)| ^ 1 and we set Tnaa{0) = Tna and Ena{0) = En 
in the following. Note that jaihri) depends on the relative displacement of the energy-level 
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compared with the two electrodes. The spin polarized current through the single junction, 
jf\H) at a=L is given from Eq.(y) with n,m ^ k'^j and G^^^{t) 9k'^ai^) 



2'7lhjso 
y<{fL{E^La)-fR{E^La)} 



(21) 



and the conductance G^\H) is expressed as G^f\H) = TLaihL)DR^(hj^) (~ 
|K-kFp-Dio-(^L)-Dija-(/i/j)). We consider the case where the left electrode is a soft magnet 
and the right one is a hard magnet. Then A^fj{h^ does not change under the inversion of 
the external magnetic field if, and the MR of the single junction is given by 



(22) 



Next we derive the MR of the double barrier structure. Conductance Gna near the n-th 
energy-level in the island is given by using df{e)/de = —6{e — jj) (T-^0) in Eq.(|l^) : 



Gna{H) 



2 • 



(23) 



[Ef - Enaihn)] + 

From this we obtain the conductance at the peak current (current at resonance) G^^^{H) 
and that at the valley current (current at off- resonance ) G°J^{H): 



GTAH)'- 
Gtm 



[rnL.(/iL) + r„^.(M]" 



27rn [EF-Ena{hnW ' 
Here we assume that Ep — Ena{hn) > 0. The change of the conductances Eq.( 
Eq.(p5l) are given in order of ^Laihi) as 



(24) 
(25) 
and 



GTAH)-GZ^'^ 

2e^ TnL^nR 



GfAH)-Gf}^ 



{YnR-TnLl^LaQlL), 



^E 



r„rr 



nL-*- nR 



^27rn [Ef-E„ 



Al,(/il)+ 



2Ena 
Ep—E„ 



(26) 



(27) 



where G'^^^^^ and GfJ-^^ are values when there is no internal magnetic field. 
Thus MR of the resonant current and off-resonant current are given as 

MRT=i: l^^^^^{ALAhL)~AU-h,)l (28) 
MRf=Y^hAUhL)-AU-hL) 

a ^ 
IE 

(7<x(/in)-7.(-/^n))]. (29) 



— Ena 

Because of the factor {TnR — TnL)/(XnR + Tni) < 1, Eq.(p8|) shows that the ratio of MR of 
the peak current is smaller than that of the single junction (Eq.(P^), whereas Eq.(pPD shows 
that the MR of the valley current is enhanced depending on the change of the energy-level 
in the island {i.e. if •jaihn) - ■ya{—hn) has the same sign as /^Laihi) - Aio-(— /il)). In the 
later numerical calculations, the enhancement of the valley current is shown as the simplest 
case where the island is non-magnetic material. Thus the MR enhancement of more than 
90% of the double junction shown by Zhang |^ is found to be the enhancement of the valley 
current and NOT that of the peak current. From these results, it is easy to conjecture the 
PV ratio, defined PVn{H) = J2aG^n^{H) G'^{H), decreases. The explicit expression of 
the PV ratio in changing the external magnetic field is given as 



PVn{H) - PVn{-H) 



a J- ni? + J- nL 

+ T7%^(7.(M-7.(-/^n))). (30) 

This shows that PV ratio is reduced when the junction conductance increases under the 
change of the direction of the external magnetic field. 

Although these results are derived starting from Eq. (pOD, they are also valid for a one 
band model where TnR^ ^ '^nRi and only one polarized spin current exists. In this case, all 
equations are similarly obtained without summation with a and the factor 1/2. 

In our formulation the resonant level is assumed to be due to a quantum dot such as 
a small magnetic particle, however, the equations derived above are general forms and the 
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results obtained here are considered to be intrinsic to the resonant tunnehng phenomenon 
such as in thin film systems ||10|| . 



III. RESULTS AND DISCUSSIONS 



Here we show the numerical results obtained from Eq.(|T8]). The current-voltage charac- 
teristics through a double barrier structure reflect the effect of DOS of the electrodes and the 



magnitude of the Fermi surface |[TT|. The form of the I-V curve has a peak when the Fermi 
energy fits the discrete energy-level and becomes lower as the bias voltage becomes higher 
in the three-dimensional electrode. The width of the peak current shows the magnitude 
of the Fermi energy when the discrete energy-level passes through the Fermi energy of the 
electrode. The existence of an internal magnetic field makes the DOS of up spins and that 
of down spins different, and the I-V curve shows a dip reflecting these two different DOS at 
the Fermi energy. These features of the model are represented in Fig. ^ which shows the I-V 
curve and MR = (J|| — J^)/J^ of a resonant tunneling structure (Eq.(p^) with that of a 
single junction (Eg. pT|)) as a function of applied bias voltage when there is a non-magnetic 
single energy-level in an island and the left electrode changes its internal magnetic field. In 
J-fl-, we take =1.8 eV, hR=2.1eV and J^, hi =-1.8 eV, hji=2.1eV for Ep = 3eV. Here we 



assume Taai^) = y ^/ EpaTaii^) for the energy dependence of the three-dimensional DOS 
of the electrodes on the I-V characteristics( Fig.l). To be more realistic, we take rQ,i(e) 
= FcoC'''^/'^^^ with constant Fqo and C . This is an extended form used in Ref. [O] and 



shows that the electron with higher energy has higher tunneling probability. As discussed 
above, the peak current is reduced more than the current of the single junction, whereas the 
off-resonant peak current is greatly enhanced, that is, more than 60%. We can obtain the 
same relation between the MR^^^ and MR^^^ even in the case where C = 1. 

The above results can also be applied to the MR of a leak tunneling current via an 
impurity trap state in the single junction [|1^]. In this case our results show that the MR of 
the leak current is less than that without a trap state and it is not necessary to take care of 
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the effect of the leak current in the measuring process of MR. 



IV. CONCLUSIONS 

In conclusion, we have studied the MR of the resonant tunneling through a double-barrier 
system by using the Keldysh formulation. Although our model is simpler than that of the S- 
matrix method by Zhang et al. , it can describe the intrinsic characteristics of the resonant 
tunneling current which the S-matrix theory cannot treat easily. We showed that the peak 
current is not enhanced by the change of the magnetization of the soft magnet whereas the 
valley current is greatly enhanced when compared with the current in a single junction. We 
also found that the PV ratio decreases when the junction conductance increases. 

T. T. is grateful to K. Sato of Toshiba Corp. for support throughout this work, and 
expresses thanks to K. Inomata, K. Mizushima, S. N. Okuno, Y. Saito, and T. Yamauchi 
for fruitful discussion. 
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FIGURES 

FIG. 1. (a) Schematic band-edge diagram for a double barrier structure. At sufficiently low 
temperature, the tunneling of electrons takes place when the discrete energy-level passes through 
the corresponding Fermi energy band, (b) Density of state of the up spin electron and down spin 
electrons at the electrode a{a = L,R). ha is an internal magnetic field. 

FIG. 2. Resonant tunneling current and the MR (see text) in the case where the left electrode 

is a soft magnet. MR^^^ is the MR of the double barrier resonant tunneling current and MR^^^ is 
that of the single barrier. In J||, Hl/Ef = 0.6 and Hr/Ef = 0.7 and in J^j, Hl/Ef = —0.6 and 
Hr/Ef = 0.7 for Ef= S.OeV where T/Ef = 10'^, rLo/EF=5.(} x 10-^ rRo/£^F=2.5 x 10"^ and 
C = 10 (sec text). The MR^^^ and MR^^^ refer to the right scale and the J|| and J|| refer to the 
left one. MR^^^ at the peak current is smaller than MR^^\ whereas that at the valley current is 
enhanced more than 60%. 
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